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This  paper  wan  prepared  by  Dr.  Shonsuke  Takagi,  Research  Physical  Scientist, 
of  the  Cold  Regions  Research  and  Engineering  Laboratory,  U.S.  Army  Terrestrial 
Sciences  Center  (USA  TSO). 

USA  ^SC  is  a  research  activity  of  the  Army  Materiel  Command. 
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ABSTRACT 

A  new  tensor  analysis,  called  the  Gibbs-Ein  stein  tensor  analysis,  is  developed  based 
oo  the  concept  that  directions  are  algebraic  quantities  subject  to  the  rule  of  forming  scalar 
products,  tensor  products,  and  linear  combinations.  The  new  tensor  analysis  is  explained 
in  this  paper  by  way  of  reformulating  continuum  mechanics  and  the  Hamiltoo-Cayley  theorem 
in  matrix  theory.  The  latter  reformulation  yields  an  explanation  of  the  deformation  dyads 
introduced  in  the  former  reformulation.  A  scalar  product  of  two  deformation  dyads  yields 
the  8 train  tensor,  which  is  a  thermodynamic  state  variable  for  thermodynamically  reversible 
deformations.  Mathematics  dealing  nth  directions  in  a  flat  space  becomes  much  simpler 
and  more  understandable  when  the  Gibbs* Einstein  tensor  expression  is  used. 


Tbl  GIBBS- HN8T1IN  TENSOR  ANALYSIS 
WITH  APPLICATION  TO  CONTINUUM  MECHANICS 
AND  CANONICAL  POEMS  OP  GENERAL  SECOND-ORDER  TENSORS 

by 

Sbunsuke  Takagi 


OfTRODOCTION 

Three  teosor  expressions  are  used  currently.  The  most  prevalent  is  the  expression  by 
eeoponeats  V1,  Vt,  TlK  T*,  etc.,  which  will  be  called  the  Einstein  expression.  The  second  expres¬ 
sion,  which  will  be  called  the  Gibbs  expression,  consists  of  linear  combinations  of  base  vectors, 
of  dyads,  of  triads,  etc.  (introduced  by  Gibbs  and  Wilson  (1901)),  whose  coefficients,  however,  are 
not  recognized  as  Einstein  expressions.  The  third  expression,  which  will  be  called  the  Gibbs- 

Einstein  expression,  is  a  combination  of  both  the  above  expressions,  expressing  a  vector  V  as 
Vt»l  m  V^,  a  second-order  tensor  Tas  TlUp j  -  T^e^  «  Tj  e,e^  «  T^eV,  a  thiid-order 

tensor  T  as  T^k»t»,»k  -  ....  etc.,  in  which  and  tf  are  covariant  and  contravariant  base  vectors, 
respectively,  defined  by 

Ve>  -  a/  (1) 

wher-  is  a  Kronecker  delta.  Coefficients  V*  and  are  contravariant  and  co variant  components 
i  •  lit  ■  »3tein  expression.  Dyads  e^,  •itf,  e^Oy  and  in  the  Gibbs  expressions  are  bases 
for  the  second-order  tensors  whose  Einstein  expressions  are  TlK  Tljt  tJ.  and  T^,  respectively. 
Similarly,  triads  ejS^...,  in  the  Gibbs  expression  are  bases  for  tbs  third-order  tensors  whose 
Einstein  expressions  are  T^k . etc. 

The  Gibbs-Einstein  tensor  expression  was  introduced  first  by  Hesseoberg  (1917)  and  ex¬ 
tended  by  tftlls  (1981).  Recently  this  notation  was  used  for  the  study  of  large  deformation  by 
Toshimura  (19577  and  Sedov  (1962)  but  it  has  not  yet  been  widely  accepted. 

Tl  1  Einstein  expression  can  bo  used  in  a  curved  space  without  introducing  normals  to  the 
curved  epace;  therefore,  it  is  convenient  for  the  study  of  iniinsic  properties  of  a  manifold.  A 
curved  sp&oo,  however,  must  be  embedded  in  a  flat  space  if  the  Gitbs-Einstein  expression  is  to  be 
applied.  This  is  because  the  differentiation  of  a  vector  belonging  to  a  curved  space  may  yield  a 
vector  that  has  as  a  conponeot  a  normal  to  the  curved  space. 

Use  of  the  Gibbs-Einstein  expression  is  based  on  the  recognition  that  identifying  directions 
with  sets  of  numbers  is  not  a  proper  definition  of  directions.  In  terms  of  axiomatic  geometry,  a 
direction  is  an  undefined  quantity,  like  a  point,  a  straight  line,  or  a  plane.  In  terms  of  abstract 
algabra,  directions  are  algebraic  quantities  subject  to  the  operations  of  scalar  product,  tensor  pro¬ 
duct,  and  linear  combinations  of  tensor  bases.  (A  tensor  product  of  vectors  is  a  juxtaposition  of 
vectors  in  a  given  order.  Vectors  in  a  tensor  product  are  non-co  mm  utative.  Juxtaposing  a  set  of 


g  '.-f-r-,  rr*  Sf'^. 
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base  vectors  forming  s  dual  basis  forms  s  sot  of  tensor  bases.  Coefficients  of  a  linear  combina¬ 
tion  of  tensor  bases, are,  in  general,  functions  of  space  and  time.)  Note  that  a  dxTereot  definition 
of  scalar  products  defines  a  different  geometry. 

The  Gibb  e-Einstein  notation  yields  simpler  expressions  end  easier  analysis  of  the  quanti¬ 
ties  containing  directions.  Geometries,  theory  of  fractions  of  many  variables,  mechanics,  and 
mathematical  physics  in  flat  spaces  should  be  reformulated  with  this  notation. 

In  tbs  first  part  of  this  paper,  the  continuum  mechanics  refonmdatad  with  the  Gibbs-Ein- 
stein  tensor  expression  will  be  summer  ii.vL  Ip  tbe  second  part,  the  Hamilton-Cayley  theorem  will 
be  reformulated  with  tbe  Gi>be-Ein  stein  tensor  etp^ssibn.  Note  that  a  matrix  to  the  Etoateia  ex¬ 
pression  of  a  second-order  tensor.  Tbe  reformulated  HamfltoCayley  theorem  to  much  simpler, 
directly  yielding  tbe  minimal  polynomial,  and  to  more  understat-dahle.  It  also  yields  a  new  con¬ 
cept  of  deforms tiju,  defining  tbe  deformation  dyad.  A  seal*  predict  of  two  deformation  dyad* 
yields  tbe  strain  tensor,  which  to  a  thermo  dynamic  state  variable  fo,  thermodynamically  reversible 
deformations. 


PARTI.  APPLICATION  TO  CONTINUIM  MECHANICS 

Let  t 1  (I  - 1,  2,  3)  be  the  coordinates  at  time  t  -  (X  A  particle  whose  initial  coordinates 
are  £*,  t*  will  be  called  particle  The  position  of  particle  tl  nt  time  t  to 


x  -  t*.  £*»  t). 

Covariant  base  vectors  mfi  «  1,  2,  3)  are  defined  by 
dx 


•i  - 


dtl 


(*) 


(*) 


Contra  variant  base  vectors  ety »  1,  2,  3)  are  defined  to  sattofy  eq  1.  Vectors  ef  represent  de¬ 
formation,  because  vector  er  for  example,  to  a  rector  obtained  by  dividing  the  rector  spanned  by 
pmtfcle*  (t1*  it K  £*,  t9)  nod  it1,  {*)  by  it1.  Jacobian  *X\  X*  X*)/*f l,  t*,  f  *>, 
where  X1,  Xs,  Xs  are  Cartesian  components  of  x.  to  equal  to  tbe  volume  of  the  parallelepfeed 


*1 x  *2  *  *r  A*8naB  tbAt  coordinates  are  right-handed,  then 


vf  -  etx  eg 


3* 


(4) 


Unit  tensor  1  to  defined  by 


i;i  which 


and 


•i  *  •/ 


g"  -  .i*  -  S*  . 


(O 


(7) 


J 


(19) 
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The  symmetric  pvt  of  gnd  v  is  denoted  by 
(gredv)*  *  «jjsV . 

Components  satisfy 


**ii  *  — —  . 

11  Dt 

Strain  tensor  «  is  given  by 


-f  ttiVVi t 
*  0 


where  ^  is  the  initial  value  of  ei  and  is  at  dependent  on  L  When  the  defamation  is  the  elonga¬ 
tion  of  «j,  Sg.  «3,  the  integral 


J* 


whose  in  teg- and  is  a  product  of  three  time-dependent  functions,  yields  a  logarithmic  strain.  In 
general,  however,  the  intepal  of  eq  23  is  dependent  on  the  path  of  integration  (Yoshimura,  1957), 
as  may  be  shown  by  following  elongations  and  rotations  in  different  aders,  and  therefae  is  not  a 
thermodynamic  state  vviable.  ■  in  eq  22  is  a  thermodynamic  state  variable  representing  a  thermo¬ 
dynamically  reversible  procoss  (see  the  end  of  this  part). 


Note  that 


a  1  it 


where  the  numbers  and  letters  unde,  the  tensor  symbols  indicate  identical  base  vectors  when  they 
are  on  different  sides  and  base  veevors  to  be  dotted  when  they  are  on  the  same  side.  The  quantity 

In  the  brackets  is  a  scalv  product  of  deformation  dyads,  and  tfi, ,  which  define  the  invase 
defamation  from  time  t  to  time  t  =  0,  To  show  this,  let  sidf<  be  a  material  point  in  tbe  neighbor¬ 
hood  of  a  particle  whose  material  bases  are  #j(  e<  at  time  t  and  V  at  time  t  =  0.  Dotting 
•jd^1  from  the  left  in  yields  e^df  Therefore,  dotting  from  the  left  in  •i9l  is  equivalent  to 
a  deformation  changing  efd^f  to  •jdf1.  In  Put  II,  the  more  realistic  interpretation  of  deformation 
dyads  will  be  given. 

Three  vectas  e^i  =  1,  2,  3)  of  f l,  f  2,  £  3  in  a  rcoce-than-three-dimensional  space  span  three- 
dimensional  subspace,  letting  f1,  £2,  £3  be  a  set  of  curvilinev  coordinates  of  the  subspace,  if 
and  only  if  and  (3tt/ are  total  differentials.  The  latter  condition,  which  yields  the 
compatibility  equations  of  components  requires  that  must  be  alinear  combination  of  vectors 
e{1  and  thaefae  shows  that  the  space  spanned  by  vectors  ef  is  Hat. 
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Let  be  a-dimensionsi  vectors,  where  nSm.  An  exterior  product  of  vectors 

(introduced  by  Gkassmann,  1844),  denoted  by  Sj  A...  A$a,  is  defined,  in  the  Gibba-Ein- 
stein  expression,  by 


a 


l 


A-1. 


(25) 


where  rr  1  is  a  permutation  symbol.  Forming  exterior  products  is  called  wedge  multiplication, 
at  wedging  for  short.  Note  that  the  rigbi-hand  side  of  eq  25  may  be  written 


»l  *« 


Jt  «  •  ••  Jl_ 
l  a 


(28) 


when  the  convention  is  applied  that  the  determinant  must  be  developed  so  that  the  elements  of  the 
first H..,mth  row  in  the  determinant  become  the  first, ....mth  base  vectors,  respectively,  in  the  ten* 
sor  products. 

The  geometric  meaning  tf  the  exterior  product  is  that 


•jA— a**  *  *v 


(27) 


where  «  is  the  tensor  expression  of  an  ^.dimensional  cube  (usually  called  orientation)  used  as  the 
unit  of  measuring  volume  V  of  the  parallelepiped  spanned  by  a} . aa. 

Expressed  by  a  three-dimensional  dual  basis  Sj,  e*  forming  a  right-handed  skew  coordinate 
system,  c  becomes 

If  V 

c  = 

-  «1/JteVe*  (28) 

where 

-  ± 

y/g 

and 

€i/a  -  y/i*uk  •  (29) 

Note  that  y/g  is  the  volume  of  the  parallelepiped  spanned  by  e^  e£,  #8. 

«  is  a  constant  tensor  fulfilling 

25  =  0  (30) 

D  t 
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f 


ft 


t 


and 


d  c 


0. 


Similarly  to  «,  lis  also  a  constant  tensor  fulfilling 


and 


D1 

Dt 


0 


dl_ 

ttl 


=  0. 


(31) 


(32) 


(33) 


An  fl-dimensional  cross  product  may  be  defined  by  dotting  with  the  D-dimensional* 


c  .  it 

123.. .0  12 


i,  io  *3 

€.  .  a  1  b  ^e  3...e  “ 

3  ...  n 


(34) 


where  ,  is  a  component  of  the  n-dimensional  « .  Because  of  the  antisymmetric  properties 

of  c ,  there  are  many  other  choices  of  dotting  base  vectors  in  c  in  the  left-hand  side  that  yield  the 
same  result  as  in  the  right-hand  side,  which,  however,  need  not  be  shown  here. 

The  exterior  differentiation  (introduced  by  Cartan  (1922))  of  a  tensor  T  of  any  order  is  given, 
in  the  Oibbs-Einstein  expression,  by 


l  n 

•  A  — 7 

di1 


(35) 


where  •i  is  usually  wedged  with  the  base  vectors  at  the  left  ends  in  the  base  tensors  of  oT/<?£\ 

but  may  be  wedged  with  any  base  vectors  in  the  base  tensors  of  dT/df  *  to  form  a  tensor  of  one 
order  higher  than  T. 

The  antisymmetric  part  of  the  three-dimensional  gradient  of  v  is  equal  to 
2  d(l 


=  —  -  *L  A 

2 \  4£l  dt*  ) 


(36) 


The  following  remark  shows  that  a  material  symmetry  that  existed  at  time  t  =  0  exists 
throughout  the  deformation. 

Remark 

Denote  t  y  e,  and  a  the  base  vectors  and  a  vector  at  time  t  =  0,  respectively, 

I  = 


(37) 
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where  scalars  a*  are  components  of  a  referred  to  e^.  At  time  t,  Sj  and  a  become  e{  and  a, 
respectively,  but  the  components  a'  are  the  same, 


a  =  »‘er 


(38) 


t  O 

Proof.  Let  initial  coordinates  different  from  f  define  fp  and  tp  different  from 
and  e, ,  respectively.  Then  we  have 


-  tpdrf 


(39) 


at  timo  t  and 


=  W 


(40) 


at  time  f «  0.  Therefore,  the  transformation  from  ei  to  fp  is  the  same  as  the  transformation  from 

h  t0  V 

Letting  a  and  a  be  one  of  fp  and  fp,  respectively,  proves  the  theorem.  The  proof  is  thus 
completed. 

The  remark  shows  that  constitutive  equations  must  be  written  in  terms  of  material 
coordinates. 

Next,  the  axiom  of  objectivity  will  be  given  the  Gibbs-Einstein  expression.  Let  cfl  »  c“ 

be  a  set  of  fixed  orthogonal  vectors  and  afl(t)  =  aa(t)  be  a  set  of  mwing  unit  orthogonal  vectors. 
Define 


Q  *  a  c°  *  a“c_. 

a  a 

12  12  12 

The  inverse  of  Q  is 


(41) 


Q-i  =  QT  *  caa“  -  c?aa 
12  12  12  ?  2 

because  they  satisfy  the  relation 


(42) 


Q.Q'l  =  Q*1  •  Q  =  1. 
la  a 2  la  a2  12 


(43) 


A  rotation  that  changes  x  =  x“ca  to  y  =  tf*ao  is  given  by 


y  =  Q  •  x  =*  x  •  Q1 


(44) 


where  nothing  is  shown  under  tensor  symbols  on  the  convention  that  two  base  vectors  adjacent  to 
the  dot,  one  on  the  left  and  one  on  the  right,  shall  be  dotted  when  no  indication  for  dotting  is 
given. 
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Adding  eq  50  and  eq  51  yields 

d'ii  +  =  q.  •i  —  +ir.  •* . qt  (52) 

d€‘  d{1  dC 

which  shows  that  (grad  v)s  is  objective.  Subtracting  eq  51  from  eq  50  shows  that  [grad  v]A  is 
not  objective. 

Base  vectors  tt  ,  e*  are  objective  as  shown  in  eq  47.  An  objective  seccnd-order  tensor 
satisfies 


T  =  tiTiitj  =  dbSbkdk 


(53) 


where  Sbk  is  the  components  in  the  moving  coordinates.  Operating  D/D t  on  T  in  eq  53  yields 


DT 

Dt 


(—  T?j  t  ip  |  \ 

\  Dt  *P  'P/  1  I 


(54) 


which  is  again  objective. 

«  in  eq  22  is  a  thermodynamic  state  variable  representing  a  thermodynamically  reversible 
process.  To  explain  this,  we  first  notice  that  d U,  for  example,  in  thermodynamics  may  be  identi¬ 
fied  with  (DI//D()Dt.  This  recognition  leads  us  to  a  thermodynamic  principle:  A  thermodynamic 

function  U,  for  example,  is  a  function  of  quantities  ql,  if  DU  is  expressed  as  a  linear  combina¬ 
tion  of  D ql  when  quantities  q*  are  independent  with  each  other. 

From  thermodynamics, 


DU  =  DQ  +  DW  (55) 

where  U  is  the  internal  energy  per  unit  mass,  and  DQ  and  DW  are  heat  and  work  inputs,  respec¬ 
tively,  per  unit  mass  per  unit  time.  Divide  DV  into  two  parts 


D V  -  (DW)rev  +  (DW)lnev  (56) 

where  (DW)rev  and  (DW)irrev  represent  reversible  and  irreversible  work,  respectively.  Then  we 
have 


DU  =  TDS  *(DW/ev 

(57) 

TDS  *  DQ  +  (DW)irr®v 

(58) 

where  S  is  the  entropy  per  unit  mass. 

When  body  couple  and  couple  stress  do  not  exist  in  the  continuum  under  consideration,  we 

have 

(DW)Tev  =  ~  o»D gir 


(59) 
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'.Therefore,  U  in  eq  57  is  a  function  of  gtj,  t,  and  S.  Because  is  a  tensor  component,  we  may 
consider  g^^  as  an  independent  variable  of  U,  where  the  Gibbs-Einstein  expression  of  eq  59 

(D! P/w  =  %/■*■  t\  ■  VJHgijVV )  (60) 

1  a  *b  b  2  12 

is  considered,,  in  which  <r  is  expressed  with  current  base  vectors.  (The  author  was  encouraged  to 
use  the  expression  on  the  right-hand  side  of  eq  60  by  Mindlin  and  Tietsten  (1962).)  «  in  eq  22  is 
integrated  to 

«  * 

-KtyW-D  (61) 

where  g^  =  ei  •  Sj  and  1  is  a  constant  tensor. 

PART  a.  CANONICAL  PORKS  OP  GENERAL  8ECOND-ORDEC  TENSORS 

The  Hamilton-Cayley  theorem  in  matrix  theory  is  given  the  Gibbo-Einstein  expression  in  the 
following.  As  shown  below,  the  dual  basis  expression  is  sore  than  suitable  for  discussing  the 
canonical  forms  of  general  second-order  tensors. 

First,  to  give  a  summary  of  this  part  and  to  show  how  the  results  may  be  used,  the  results 
will  be  applied  to  three-dimensional  tensors.  Canonical  forms  of  not  necessarily  symmetric  real 
three-dimensional  second-order  tensors  in  the  Gibbs-Einstein  expression  are  classified  into  four 
Categories: 

Category  1:  Eigenvalues  A,(i  =  1,  2,  3)  are  all  real,  and  determine  three  pairs  of  left  and  right 
eigenvectors  which  are  never  orthogonal  with  each  other.  Then  eigenvectors  span  a  dual  basis 
S|,  e1  satisfying 

e1  •  T  *  Aj  •l  (62) 

T  •  •*  -  \ttl  (63) 

where  i  =  1,  2,  3.  The  summation  convention  is  not  jpplied  on  the  right-hand  sides.  Juxtaposing 

e(  from  thB  left  in  eq  62  and  e{  from  the  right  in  eq  63  with  tne  summation  convention  applied 
yields  the  same  expression 

3 

T  -  £  Af  tl*t  (64) 

which  is  the  canonical  form  of  category  1.  Note  that  eigenvectors  are  not  necessarily  unit  nor 
orthogonal.  When  e,  *  e*  eigenvectors  are  unit  orthogonal  and  T  is  symmetric.  Eigenvalues  in 
category  1  may  be  multiple  roots. 
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Category  2:  A  ^  is  a  double  real  root  which  determines  an  orthogonal  pair.  An  orthogonal  pair  is 
a  pair  of  left  and  right  eigenvectors  determined  for  an  eigenroot  and  orthogonal  with  each  other. 
Choose  e2  and  e1  as  the  left  and  right  eigenvectors  determined  for  Ajt  then 

e2  '  T  =  AjOg 

(65) 

T  •  e1  =  Aj  e1. 

(66) 

Let  A2  be  the  remaining  real  eigenroot,  and  e3  and  e3  be  the  left  and  right  eigenvectors,  respec 
tively;  then, 

e3  •  T  .  A2e3 

(67) 

T  •  e3  =  Ag*3. 

(68) 

As  proved  later,  and  e2  can  be  chosen  to  satisfy 

•  T  =  A(ej  +  «2) 

(69) 

and 

T  •  e2  =  A(e1  +  e2). 

(70) 

Juxtaposing  e1,  e2,  e3  from  the  left  in  eq  69.  65,  67,  respectively,  and  summing  the  results,  and 
juxtaposing  e.,  e2,  e3  from  the  right  into  eq66,  70,  68,  respectively,  and  summing  the  results, 

yield  the  same  expression. 

T  =  A^e^j  +  »2«2  +  e^g)  +  A2ese3 

(71) 

which  is  the  canonical  form  of  category  2.  The  canonical  form  has  one  off-diagonal  term  AjO^Sg. 

As  shown  later,  choice  of  base  vectors  for  expressing  the  canonical  form  in  category  2  is  not 
unique. 

Category  3:  A  is  a  triple  real  root  which  determines  an  orthogonal  pair.  Chose  e3  and  e1 
as  the  left  and  right  eigenvectors  determined  for  A;  then 

e3  •  T  =  Ae3 

(72) 

• 

>-* 

ii 

fcr 

• 

(73) 

As  proved  later,  elt  e2.  e2  and  e3  can  be  chosen  to  satisfy 

0^1  =  A(el  +  e2) 

(74) 

V  T  =  ^e2  +  e3) 


(75) 
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T  •  s2  =  Ate1*#2) 


T-s3  =  A(e2  +  e3) .  (77) 

Juxtaposing  e1,  e2  e3  from  the  left  ia  eq  74.  75.  72,  respectively,  and  summing  the  re¬ 
sults.  and  juxtaposing  er  «2>  «3,  from  the  right  in  eq  73.  76.  and  77.  respectively,  and  summing 

the  results,  yield  the  same  expression, 

T  =  +  e*e2  +  e2e3)  (78) 

which  is  the  canonical  form  of  category  3.  The  canonical  form  has  two  off-diagonal  terms,  Ae^g 
and  Ae2«3.  As  shown  later,  choice  of  base  vectors  for  expressing  canonical  forms  in  category  3 
is  not  unique. 

Category  4:  Two  eigenroots  are  conjugate  complex.  As  shown  later,  eigenvalues  and  eigenvec¬ 
tors  in  this  case  are  expressed  as  p( cos  d  *  >  sind)  and  e}  t  i«2,  e1 1  fe2;  thus  we  find 

(«j  ♦  i«2)  •  T  =  p(cos 0  *  isindKSj  t  i«2)  (79) 

T  •  (e1  ie2)  *  p( cos  d  ?  /  sindMe1  v  f  e2)  (80) 

where  «2>  s1,  s2  are  real  vectors,  and  p  and  d  are  real  numbers.  Deconposing  eq  79  and  80 
into  the  real  and  imaginary  parts  yields 

•j  •  T  »  p(*j  oos  d  +  *2  sin d)  (81) 

•  T  *  p(®2  ccs  d  —  e jSin  d)  (82) 

T  •  a1  *  pfe1  cos  d  -  e2 sind)  (83) 

T  •  s2  *  p(e2cos  d  +  «l  sind) .  (84) 

Let  ft  be  the  remaining  real  eigenroot  and  e3,  e3  be  the  left  and  right  eigenvectors,  respectively; 
then, 

•3  •  T  *  P®3  (85) 


T  •  s3  =  ftt3  .  (86) 

Juxtaposing  e1,  e2.  e3  from  the  left  in  eq  81.  82,  85.  respectively,  and  summing  the  results,  and 
juxtaposing  Sj.  «2,  e3  from  the  right  in  eq  83  ,  84,  and  86,  respectively,  and  summing  the  results. 

yield  the  same  expression 
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T  -  p^e^j-i- e2e2)cos0  +  p(n1n2-n2e1)sin0  +  pn3n8  (87) 

which  is  the  real  canonical  form  of  category  4. 

A  general  second-order  tensor  T  can  be  written  as  a  dyad  T  «  a1^  and  may  be  interpreted 
as  a  deformation  dyad.  Interpreted  this  way.  the  equations  in  eq  62  show  that  elongations  in  three 
directions  er  e^  and  e3  have  occured.  Equation  6?  in  category  2  shows  that  a  slip  has  occurred 

in  the  a1(  eg  plane  along  the  «2  axis.  Equations  74  and  75  in  category  3  show  that  a  double  slip 
has  occured  in  the  e1#  #2  plane  along  the  eg  axis  and  in  the  #2,  e3  plane  along  the  a^azis. 

Equations  SI  and  82  in  category  4  show  that  a  rotation  by  angle  0  has  occured  with  a3  as  the 
axis  of  rotation. 

Complex  tensors 

In  the  following  sections,  general  e-dime neional  second-order  tensors  are  given  canonical 
forms.  For  that  we  must  first  extend  theory  of  real  tensors  to  theory  of  complex  tensors. 

A  set  of  unit  orthogonal  vectors  cq  -  c“  (a  -  1.....&)  is  fixed  in  the  space  aM  used  as  the 
standard  of  the  coordinate  systems.  Vector 

v  -  v“ca  (88) 

is  called  a  complex  vector  if  conponents.  v°  *  va  («  -  l,...,a)  in  the  standard  expression  (eq  88) 
are  complex  numbers.  Conjugate  v  of  v  is  defined  by 

v  -  v°ca .  (89) 

Dotting  O  complex  vectors  n  with  v.  denoted  by  n  -  v,  is  defined  by 

n*F  -  n  -  v  (90) 

where  dotting  ( ♦)  on  the  right-hand  side  is  the  dotting  in  real  Euclidean  geometry  applied  to  complex 
vectors.  Vector  v  satisfying  v  -  v  *  1  is  said  to  be  of  unit  length.  Vectors  n  and  v  satisfying 
■  -  v  -  0  or  v  •  n  -  0  are  said  lo  be  orthogonal.  A  dual  basis.t{.  tf  for  complex  vectors  satisfies 

V  #  -  sJ  (9» 

or 

«  8*1  («) 

where  and  SJj  are  Kronedcer  deltas. 

The  standard  expression  of  an  r- dimensional  complex  secood-arder  tensor  T  is  defined  by 
T  -  T/n-'n, 


(93) 
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i1  =  ifi'i  (105) 

Proof.  The  foregoing  two  lemmas  are  proved  in  the  following 

Dotting  e1*,  or  efe'  from  the  right  in  victor  v  expressed  >,s  or  v ^  proves  that  1  for 
the  right  multiplication  is  is^  or  e^,  respectively.  Dotting  e'e^  or  ip'  from  the  left  in  vector  v 
expressed  as  v^V  or  improves  that  1  for  the  left  multiplication  is  also  or  e^e1.  The  proof 
that  e'Oj.  =  tp1  follows. 

Juxtaposing  t?  from  the  right  in  eq97  yields  eq  102.  Equations  103,  104  and  105  are  de¬ 
rived  similarly.  Substituting  eq  105  intoe^ey  yields 

•S  =  (a) 

Substituting  eq  103  into  e^1  yields 

Bj*'  =  g'^Vy .  (b) 

Comparing  eq  a  and  b  proves  a  part  of  eq  99.  The  rest  of  eq  99  is  proved  similarly.  The  proof  is 
thus  completed. 

Define  the  transformation  of  base  vectors  by 

=  ataca  (106) 

and 

ef  =  fPbJ  (107) 

where  aj  and  bj  are  complex  numbers  satisfying 

a fbj[  =  (108) 

Equations  106  and  107  are  transformed  in  the  following  to  their  inverses. 

Lemma  lc 

c“  =  e'a/1  (109) 

and 
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Proot.  Dotting  eft  into  eq  106  yields 

aft  =  zl  •  eft  .  (a) 

Juxtaposing  from  the  left  of  eq  a  yields  eq  109.  Equation  107  transforms  to  eq  110  similarly. 

Theorem  1 


«iabJ~c, 

(111) 

Va  “b‘ . 

l  u 

(112) 

Proot.  Substituting  eq  110  or  109  into  eq  ICS  or  107  yields  eq  111  or  112,  respectively. 

Corollary  1 


-  ifiaS  (113) 

•i  ■  A6  V  •  (li4) 

Proof.  Taking  the  conjugates  of  eq  111  or  112  proves  eq  113  or  114,  respectively.  The 
proof  is  thus  completed. 

Equations.  Ill  through  114  ate  the  equations  that  must  be  used  to  transform  non-standard 
expressions  to  standard  expressions. 

Eigenvectors  and  eigenvalues 

Let  T  be  an  n-dimensionai  tensor  in  the  standard  form.  A  left  or  right  eigenvector  of  T 
is  an  n-dimensional  complex  vector  x  or  y  such  that  dotting  ( • )  x  or  y  from  the  left  or  right  in  T 
yields  a  vector  in  the  direction  of  x  or  y,  respectively;  that  is 

x  •  T  -  Ax  (115) 

or 

T  y  =  Ay  (116) 

where  A  is  an  eigenvalue. 

Let  T  and  a  left  eigenvector  x  do  expressed  as 

T  •=  T/e'e;  (117) 

and 

x  -  x\  (118) 

respectively.  Substituting  eq  117  and  118  into  eq  ;>5  and  equating  the  components  on  both  sides 
of  the  transformed  equation  yields 
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x1  T/  =  Xx’ .  (119) 

Let  a  right  eigenvector  y  be  expressed  as 

y  =  ybtb.  (120) 

Substituting  eq  117  and  120  into  eq  llo  yields 

T/  9j  =  Ay,..  (121) 

Equations  119  and  121  determine  the  same  characteristic  equation 

|  TJ  -  AS/  I  =  0.  (122) 


Therefore,  o*-e  A  determines  at  least  one  left  and  one  right  eigenvector.  The  rank  of  matrix 
(T/  -  AS/)  will  be  Celled  the  rank  of  root  A.  Solutions  of  eq  115  and  116  do  not  depend  on  the 
choice  of  a  type  of  T  or  of  a  dual  basis. 

Canonical  forma 
Lemma  2a 

Let  A  and  n  be  two  different  eigenvectors  of  T;  let  x  and  y  be  left  and  right  eigenvectors 
determined  for  A,  respectively;  and  let  u  and  v  be  left  and  right  eigenvectors  determined  for  /r,  respec¬ 
tively.  Then  x  and  v  are  orthogonal. 

x-v  =  0  (123) 

and  y  and  a  are  orthogonal, 

y  •  n  -  0  .  (124) 

Proof.  Dotting  (• )  v  from  the  right  in  eq  115  yields 

x  •  T  •  v  «  Ax  •  v.  (a) 

Dotting  (• )  from  the  left  in  the  equation 
T  •  v  =  fiv 

yields 

x  •  T  •  v  *  /tx  •  v.  (b) 

Because  A+  n  by  assumption,  eq  a  and  b  are  compatible  if  and  only  if  eq  123  is  true.  Equation  124 
may  be  proved  similarly. 

Lemma  2b 

If  all  the  roots  are  single,  no  orthogonal  pair  exists. 


!?7T  IVKWBW 


ffm  u  w 
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The  next  lemma  proves  this  lemma.  Note  that  vectors  ef  and  tl  forming  a  dual  basis  are 
not  orthogonal. 

Lemma  2c 

Assume  that  n  roots  Aj . An  of  the  characteristic  eq  122  are  distinct.  Then,  T  can  be 

reduced  to  the  canonical  form 

a 

T  .  y  A (125) 

tu-i  1  1 

i-1 

The  rank  of  A t{l  =  1,..., n)  is  equal  to  n  -  1. 

Prool.  Each  Ai  determines  at  least  one  left  eigenvector.  Let  be  the  left  eigenvector 
determined  for  Ai  (i  ■  l,...,n) 

VT-Vt  <a> 


V  T  *  Va*  <b> 

Form  a  dual  basis  ef,  ei  (i  =  l,...,n).  Juxtaposing  e* . e°  from  the  left  in  eq  a,...,  eq  b, 

respectively,  and  summing  the  result  yield  eq  125.  Dotting  (• )  tl  from  the  right  in  eq  125  shows 
that  e1  is  a  right  eigenvector  determined  for  A {(i  »  1, 

T  in  eq  125  shows  that  each  Af  determines  one  and  only  one  pair  of  left  and  right  eigen¬ 
vectors.  The  rank  of  the  roots  is  therefore  all  equal  to  n  -  1.  The  proof  is  thus  completed. 

A  multiple  root  which  determines  fewer  eigenvectors  than  the  number  of  multiplicity  of  the 
root  is  said  to  be  singular. 

Lemma  2d 

At  least  one  orthogonal  pair  belongs  to  a  singular  multiple  root.  If  more  than  one  orthogonal 
pair  belongs  to  a  singular  multiple  root,  the  number  of  orthogonal  pairs  can  be  reduced  to  one. 

Proof.  Assume  that  eigenvectors  a  and  y  determined  for  a  multiple  root  are  not  orthogonal. 
Then,  we  can  choose  ev  *  x  and  e1  =  y  after,  if  necessary,  changing  the  lengths  of  x  and  y  to 
make  x  *  y  »  1 .  Then  we  have 

•  T  *  ACj  (a) 

and 

T  •  e1  *  Ae1 .  (b) 

Forma  dual  basis  •l,  e<  (/  -  l,...,n)  by  introducing  a  certain  number  of  base  vectors  that  are 
independent  of  e^,  eland  each  other;  express  T  with  the  dual  basis  'ef,e‘ 
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T  =  TjPtj.  (c) 

Substituting  eq  c  into  eq  a  and  b  yields 
Tj1  =  A 

T*  =  ...  =  T°  =  0 
Tg1  =  ...  =Tn1  =  0. 

Therefore,  T  in  eq  c  becomes 

T  =  Ae^+S  (d) 

where 

S  =  TpVe,  (e) 

in  which  p  and  q  represent  2 . n.  8  is  an  (n  -  l)-dimensional  tensor  such  that  one  of  the  eigen¬ 

values  is  A.  Therefore,  if  the  number  of  non-orthogonal  pairs  is  equal  to  the  number  of  multiplicity, 
eigenvectors  determined  for  the  multiple  root  can  form  a  dual  basis,  which  shows  that  at  least  one 
orthogonal  pair  must  belong  to  a  singular  multiple  root. 

If  more  than  one  orthogonal  pair  belongs  to  A,  choose  two  pairs  x,  y  and  n,  v.  Form  two 
non-orthogonal  pairs  by  setting 

=  x,  e1  =  v 

e2  =  u,  e2  =  y. 

Then,  T  can  be  reduced  to 

T  =  Ae^j  +  Ae2eg  +  R  (f) 

where 

R  =  Trseres  ® 

in  which  r  and  s  represent  3 . n.  Continuing  this  process  the  number  of  orthogonal  pairs  belong¬ 

ing  to  A  can  be  reduced  by  an  even  number,  but  cannot  be  reduced  to  less  than  one.  The  proof  is 
thus  completed. 

Theorem  2 

Assume  that  A  is  a  (t  +  r)-tuple  root  of  rank  n  -  r  -  1  of  an  n-dimensionai  tensor  T.  Then, 

T  can  be  reduced  to 


UK  t  JH  J-J  U'J’TVJrT 
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where  8  does  not  have  X  as  an  eigenvalue  and  is  (n  -  t  -  ^dimensional.  Then,  the  left  eige> 
vectors  determined  for  X  are  •t+1»....«t+l.  The  right  eigenvectors  determined  for  X  are 

el,  et+l,...,et+r.  Vectors  ^ . et  satisfy 

VT  =  AV  +  Vi  (127) 

where  <f>  represents  t  -  1  integers 
2  <$<t. 

Vectors  e1,...,®1”1  satisfy 

T  •  ®^  -  Xe^  +  «^+ 1  (128) 

where  <f/  represents  t  -  1  integers 
1  <  ^  <  t-  1  . 

The  number  of  off-diagonal  terms  is  t  -  1. 

Proof.  Because  the  rank  of  X  is  n  -  r  -  1,  (r+  1)  pairs  of  independent  left  and  right  eigenvectors 
exist,  of  which,  if  tS  2,  one  is  an  orthogonal  pair.  Denote  by  ®r  tl  the  orthogonal  pair,  by 
Oj,  e{+1 . et+r  the  left  eigenvectors,  and  by  e£,  et+1 . e£+r  the  right  eigenvectors.  Then 


®i  ‘  T  * 

x«i 

(a) 

et+rT 

=  Aec+i 

(b) 

®t+r  *  T 

=  Aem 

(c) 

T  •  ®£  = 

X®£ 

(d) 

T  •  ®£+1 

=  X®£+* 

(e) 

T  •  ®£+r 

• 

=  Xet+r. 

(0 

Form  a  dual  basis  «f,  «£ (i  =  l,...,n)  by  arbitrarily  introducing  independent  vectors  ®2,...,  «t, 
•t+r+f  •*■•»■  •1-— *t-1>  «t+r  1 . Express  T  as 

T  =  T/iiej  (g) 
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Subsiding  eq  g  into  eq  a  through  f  yields 


Ti  - 

■  ^ 

(a1) 

w 

= 

(b') 

Tj 

= 

(c1) 

Tt 1  = 

AS/ 

(<!') 

T/+1 

=  AS/+1 

(e'J. 

T/+f 

=  AS/+'  . 

(f) 

The  matrix  of  T  -  A1  is  shown  in  Figure  1,  in  which  all  the  elements  on  the  straight  lines 

are  zero,  submatrix  A  is  composed  of  2nd . tth  row  and  1st . (t-  l)th  column,  and  the  main 

diagonal  of  A  is  not  on  the  main  diagonal  of  the  matrix  of  T  -  Al. 

Define  x  fa  a  given  vector  a  by 

x  •  T  =  Ax  +  a.  (h) 


Figure  1.  Matrix  of  T-Al. 
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Equation  h  becomes 

x*(T/-AS/)  -  a1  (i) 

by  substituting  eq  g  for  T,  eq  118  for  x,  and  a  *  . 

Because  of  conditions  in  eq  d'  through  f\  equations  in  eq  i  are  simultaneous  n  -  r  -  1 
equations  given  for /=.  1,...,  t-1,  t+r+l,...,n.  Vector  a  must  satisfy  a  condition  that  components 

a1,...,  at+r  are  zero.  Because  of  conditions  in  eq  a*  through  c\  equations  in  eq  i  have  n  -  r  -  1 

unknowns  x2 . xl,  xt+r+1,...,  xn  .  The  determinant  of  the  simultaneous  equations  (eq  i)  is  of 

rank  n  -  r  -  1,  and  is  the  only  non-zero  (n  -  r  -  1>  dimensional  submatrix  containing  A.  Then  x 

is  determined  with  (r  +  l)  arbitrary  components  x  *,  xt+  xt+r. 

Put  a  «  and  define  x2  by 

x2  T  =  Ax2  +  er  (j) 

xg  is  linearly  independent  of  e1.  To  show  this,  define 


y  »  yl»l  +  yzxz . 


00 


We  find 

y-T-Ay  +  y2*!  (1) 

by  use  of  eq  a  and  j.  Therefore,  if  y  -  0,  we  necessarily  have  y1  =  y2  =  0.  Let  xg  be  ctosen  as 
new  eg.  Then  eq  j  becomes  one  of  eq  127. 

Put  a  =  e2  and  define  x3  by 

*3‘T~a*8  +  V  (">) 

x3  is  linearly  independent  of  e2  and  «lt  as  may  be  shown  similarly  to  the  above.  Let  x3  be 
chosen  as  new  e3 .  Then  eq  m  becomes  one  of  eq  127. 

Continuing  this  process,  vectors  e4,...t  et  are  defined  and  eq  127  is  proved.  All  the  vectors 
represented  by  e^_j  in  eq  127  satisfy  the  condition  which  must  be  satisfied  by  vector  a  in  eq  h. 

Suppose  that  eq  g  is  expressed  with  the  new  base  vectors  thus  introduced.  Substituting 
eq  b  thus  determined  into  eq  127  shows  that  matrix  A  in  Figure  1  is  a  unit  matrix  and  that  sub¬ 
matrix  C  in  Figure  1  is  a  zero  matrix,  yielding  eq  126. 

Dotting  ( • )  from  the  right  into  T  in  eq  126  yields  eq  128. 

Substituting  eq  g  into  eq  128  shows  that  submatrix  D  in  Figure  1  is  a  zero  matrix.  The 
proof  is  thus  completed. 
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Integer  t  is  called  the  grade  of  A  (Turnbull  and  Aitken,  1932).  Vectors  «2,.. 

et  are  called 

left  pseudo-eigenvectors  determined  for  A.  Vectors  e1,...,  et_1  are  called  right  pseudo-eigenvectors 
determined  for  A.  T  -  S  in  eq  126  is  called  the  canonical  part  belonging  to  A. 

Corollary  2a 

The  canonical  part  belonging  to  A  does  not  change  form  when  dual  basis  e(, 
formed  to  another  dual  basis  defined  by 

e*  is  trans- 

-  P,ei 

(129) 

f2  =  p\*P% 

(130) 

f£-1  =  pt”1e1  +  pt"2e2+...  +  P  xet_  t 

(131) 

t  t+r 

ft  =  ^  Pt+1"S  +  ^  ^ 

(132) 

t+r 

*t+i  =  2 

£=t+i 

(133) 

t+r 

*t+r  =  tt+rel+  ^  W‘4 

f=t+i 

(134) 

where  p1.....?*,  if,  fj,  tj  U  <  =  t  +  1 . t+r)  are  arbitrary  if  they  make  fr. 

independent. 

...ft+r  linearly 

Proof.  If  vectors  f  j . f(+f  satisfy 

fj  •  T  =  Afj 

(a) 

VT  =  AVf*-i’  (2  =  *=<t) 

(b) 

f^-T  =  Al^,  (t  +  U  t+r) 

(c) 

the  canonical  part  determined  by  lj . ft+r  ia  the  same  as  that  determined  by 

Let  index  a  represent  integers  l,...,t+r.  Substituting 


and  T  in  eq  126  into  eq  a  shows  that  vector  f ,  must  be  such  as  in  eq  129. 
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Substituting 


and  T  in  eq  126  into  eq  b  yields  the  conditions, 

,  2  _  /  1 

%  3  Vi 


/  t  _  /  ‘-1 

V  -  V-1 


Vi* 


Vi+r  = 0 


which  snow  that  vectors  f„ . f£  must  be  such  as  in  eq  130 . eq  132,  respectively-  Substituting 


and  T  ineq  126  into  eqc  shows  that  vectors  f£+1 . flrl  must  be  3uch  as  in  eq  133 . eq  134. 

The  proof  is  thus  completed. 

Corollary  2b 

Left  aud  right  eigenvectors  determined  for  a  single  toot  are  not  orthogonal. 

Pmol.  Denote  by  8  the  sum  of  all  the  canonical  parts  belonging  to  the  multiple  roots.  Uten, 
all  the  eigenvalues  of  T  -  8  are  single,  and,  as  proved  by  Lemma  2b,  no  orthogonal  pair  exists  for 

T-8. 


Haadlt'Mfr-Cayiey  Theorem 


Let  Oj . ea  be  n  independent  vectors  and  8  be  an  n-dimensiocal  complex  tensor  of  any 

order.  If  the  relation 

e,  •  8  =  0  (135) 

is  true  for  all  e^i  =  l,...,n),  then  8  is  identically  equal  to  sera 

Proof.  Form  a  dual  basis  e;,  e‘.  Juxtaposing  s'  from  the  left  into  eq  135  and  summing 
over  i  yields  the  required  property  of  8.  The  proof  is  thuscroplcted. 

Let  T  be  a  complex  second-order  tensor  in  a  standard  form 

T  =  e1  tj  .  (136) 

Then,  a  time  scalar  product  of  T 


T 


T  =  T* 


(137) 
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abbreviated  to  T*.  is  a  second-order  tensor 


-ijji  j  i» 

i  ix 


where  a  is  any  positive  integer. 


(138) 


Let  \u  (1/  =  be  the  roots  of  the  characteristic  equation  of  T.  Assume  that  A^  is 

(tu  +  r^Mqple  root  of  rank  a  -  rv  -  1,  where 

p 

^  (t„  ♦  rv)  =  s  .  (139) 

v  =  1 

Then  a  second-order  tensor 

f(T)  =  n  (140) 

v=  1 

is  identically  equal  to  sera 

Proof.  Let  tbe  collection  of  eigenvectors  and  pee  jdo-fcigen.ee tors  belonging  to  all  the 
eigenvalues  form  a  dual  basis  ef 


Equations  for  A{  are 

•j  •  (T -  Ajl)  *  0  (a) 

Sg-tT-Ajl)  =  (b) 

•  • 

•cJ'(T-Ai>)  =  •tj-l  <c> 

eti+I  (T-Ajl)  =  0  (d) 

•t^r,  •  (T-A.l)  -  0.  (e) 

Equation  b  through  c  can  be  changed  to 

Sg  •  (T  -  Ajl)2  =  0  (0 

e.  (T-A  I)*1  =  0.  (g) 

ll 
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Therefore,  we  find 


{ 


(T-a  l)*1  =  0 


where 


Similarly.  we  find 

(T-A  i)*"  =  0 


(b) 


(i) 


^  J  1 . 2  vv- 

/t=l  |t=l 

Fora  A  T)  ha  eq  122.  The  order  of  dotting  in  eq  122  any  be  exchange*:,  because  the  re¬ 
acts  of  dotting  in  differa*  orders  are  the  line.  Vie  therefore  find 

•,  -  N't)  «  0  (j) 

where  i  «  .  Then.  AT)  awat  be  identically  eqaal  to  aero.  The  proof  is  tbns  coapleted. 

Polynomial  A*)  obtained  by  submitting  scalar  x  far  T  is  the  ainieal  polynomial  defined 

for  T. 


A  tenor  is  said  to  be  teal  if  dare  exists  a  transformation  that  changes  all  the  conponeics 
to  real  nnaters  and  all  the  bane  vectors  to  real  vectors  at  the  saae  tine.  In  the  following.  T  is 
a  second-order  tensor  that  has  real  compoeects  and  real  base  vectors. 

It  is  obvioas  that  a  real  ejgenvabe  of  T  determines  real  eigenvectors,  and  that  coejngate 
conptexeigemralnse  of  T  determine  coojngate  conplex  eigenvectors.  Only  the  latter  case  need 
be  discus r»d. 

As  some  tfcx  a  tifi  are  (t  +  r)-ta|4e  roots  of  rank  n  -  r  -  1  of  a  real  secood- order  tensor  T. 

she  re 


2lt  «-  r)  5  a  . 


Then.  T  can  be  reduced  to 
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2  (t+r)  t+r 

T  =  a  £  fPfp  +  J8  2(f2,f29-l-f2,‘lf29)  + 

P-1  <P=.i 

i 

+  f80‘If2*-3)+S  U«> 

where  f  ,  P-’  [  1  £  p  |  2(t  +  r)]  are  real  vectors,  forming  a  dual  basis,  and  8  does  not  have  a  tip 


as  eigenvalues  and  is  2(n  - 1  -  r>dimensional. 

Vectors  tp,  tp  (p  =  l,...,2(t+r))  satisfy 

f1-T  =  af1-^f2  (142) 

fg  •  T  =  af2  +  /3f  1  (143) 

*20-1  '  T  =  a*20-l  “  ^*20  +  *20-3  (144^ 

f20  ‘  T  =  a*20  +  /3*2 0- 1  +  *20-2  (2  £  0  £  t)  (145) 

*2£-rT  -  a*2«f-l-0*2<f  <146> 

*2^ '  T  =  a*2f +  ^*2f  l  (l  + 1  <  f  ^  t +r)  (147) 

T  •  1  =  at^Y~ 1  +  /Sf2>'  +  f2^ 1  (148) 

T  •  &  =  (l<y^t-l)  (148) 

T  •  f27-l  =  aW-l  +  pi*]  (150) 

T  •  f27?  =  atZri-ptZrl'x  (t$7]<t+r).  (151) 

Complex  vectors  —  U2q.l  ±  it2q) .  ~U2q* 1  +  if24)  (5  =  1 . t+r)  form  a  dual  basis 

s/2  V2 

Proof.  Let  Cj . e{4f  be  complex  eigenvectors  and  pseudo-eigenvectors  belonging  to 

a  +  i/3 

«J  •  T  =  (a  +  i/3)e  x  (a) 
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eg  •  T  =  (a  +  i/3)e  x 


e£  •  T  =  (a  +  i/9)et  +  «t_  j 
®t+i  '  T  =  (a  +  i£)«t+1 


®t+r  •  T  -  («  +  ^>®r+r  •  (e) 

Write  efl(l  ^  g  I  t  +  r)  with  real  vectors  f2?  as 

®g  =  -^jr  ^aq-i +  ^2q)'  (0 

Substituting  eq  f  into  eq  a  through  e  yields  eq  142  through  147. 

Determine  fp[  1  £  p  ^  2(t+r)]  to  form  a  dual  basis  tpl  fp.  Juxtaposiag  f1 . f2(t+r)  from  the 

left  into  eq  142  through  147.  respectively,  and  adding  the  equations  thus  formed  yield 

T  =  f1taf1-/3f2)  +  f2(af2  +  ^f1)  + 

t 

+  p2<?i-1{af2<£-i~/3*2<£  +  hf-a)  +  f2<^^“?20  +  ^f20-l  +  *2c£-2^1  + 

<f>=2 

t+r 

+  2  [(2^' 1  (flf2^-  r  ^  +  ^ (at2f +Pl2£- 1>1  +  s  (g) 

f-‘+l 

where  S  does  not  have  a  tip  as  eigenvalues.  T  in  eq  g  oau  be  transformed  to  T  in  eq  141. 

Dotting  t L,...,f2(t+r)  from  the  right  in  T  in  eq  g  yields  eq  148  through  151,  respectively. 

The  proof  is  thus  completed. 

T  in  eq  g  may  be  expressed  as 


T  =  8+  [fJ...  f2<t+r>]T  .' 


2(t  +  r) 


where  T  is  a  2 (t+r)  by  2(t+r)  matrix 


D  BY 


in  which  A  is  a  2t  by  2t  matrix 
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A 


a  -0  0  00 . 000  0 

0  a  000 . 000  0 

1  0  a-aO . 000  0 

Ol0aO . 0000 

0  0  1  0  a . 0000 


0  0  0 

0  o  •  •  • 

•  •  I  0  -«  -0 

0  0  0 

o 

O 

••010“ 

B  is  a  2r  by  2r  matrix 


i 

i 

i 

i 


'a  -0  0 
0  a  0 
0  0a 


0  0  0 
0  0  0 


oo  o: 
0  0  0 
0  0  0 


O  a  -0 
0  0  a. 


C  and  0  are  composed  of  zero  only.  In  matrix  A,  all  the  elements  on  the  line  parallel  to  the  main 
diagonal  and  containing  1  are  1. 
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